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ABSTRACT 

We present velocity dispersion measurements of 20 Local Group stellar clusters (7 < log(age [yrs]) < 
10.2) from integrated light spectra and examine the evolution of the stellar mass-to-light ratio, T*. We 
find that the clusters deviate from the evolutionary tracks corresponding to simple stellar populations 
drawn from standard stellar initial mass functions (IMFs). The nature of this failure, in which is at 
first underestimated and then overestimated with age, invalidates potential simple solutions involving 
a rescaling of either the measured masses or modeled luminosities. A range of possible shortcomings 
in the straightforward interpretation of the data, including subtleties arising from cluster dynamical 
evolution on the present day stellar mass functions and from stellar binarity on the measured velocity 
dispersions, do not materially affect this conclusion given the current understanding of those effects. 
Independent of further conjectures regarding the origin of this problem, this result highlights a basic 
failing of our understanding of the integrated stellar populations of these systems. We propose the 
existence of two distinct initial mass functions, one primarily, but not exclusively, valid for older, metal 
poor clusters and the other for primarily, but not exclusively, younger, metal rich clusters. The young 
(log(age [yrs])<9.5) clusters are well-described by a bottom-heavy IMF, such as a Salpeter IMF, while 
the older clusters are better described by a top-heavy IMF, such as a light-weighted Kroupa IMF, 
although neither of these specific forms is a unique solution. The sample is small, with the findings 
currently depending on the results for four key clusters, but doubling the sample is within reach. 
Subject headings: stars: formation, luminosity function, mass function; galaxies: fundamental param- 
eters, evolution 



1. INTRODUCTION 

The stellar mass of a galaxy is intricately connecte d 
to that galaxy's e nvironment (|K auffmai m et al.l 120041 ). 
metal abun dance (iTremonti et all [2004), star forma- 
tion history dBundv et aL| | 2006D . dark matter halo mass 
(|van den Bosch et al.1 120031) . and just about any other 
principal characteristic one cares to consider. Yet, we 
rely on estimates of stellar masses that are potentially 
rife with systematic errors. Those estimates hinge on 
our poor knowledge of a galaxy's star formation history 
and are further predicated on two aspects of stellar evo- 
lution that are poorly understood: stars' behavior during 
the phase(s) of their life at which they are at their most 
luminous and the initial distribution of stellar masses 
(the initial mass function or IMF) . Despite these various 
complications, we utilize simple models to estimate stel- 
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lar masses because dynamical measurements are techni- 
cally difficult to obtain and complicated by the presence 
of dark matter. Even when a dynamical mass estimate 
is available, modeling is required to determine the rela- 
tive contributions of stars and dark matter to the total 
mass. Given that any apparent failure of such models 
impugns our understanding of either stellar evolution or 
the IMF, and therefore affects many aspects of our study 
of the extragalactic universe, testing these models on the 
simplest possible stellar systems is paramount. 

Discrepancies between models and observations, when 
found, are commonly attributed to deviations in the 
IMF from th e adopted prescrip t ion (for some recent 
exam p les see | van Dokkuml 120081: iDabringhausen et al.l 



2001T iTreu et al.l 120101: Ivan Dokkum fc Conrovl 120101: 
Dutton et al. 1 120121 ) rather than to either the star for- 
mation history or stellar evolutionary models. Direct 
measurements of the ini tial mass function a re difficult 
for various reasons (see iBastian et al.l [2010, for a re- 
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view), particularly over the full range of environments 
and conditions. Various "standard" descriptions of the 
IMF are often used to span the range of possibility al- 
though given the simple mathematical descriptions used 
for the IMF it is not evident that such an approach is 
complete. Some of the IMFs in widest use incl ude the 
original power-law description by iSalpeterl ([1955D , which 
is generally acknowledged to be a poor fit to current 
data f or nearby regions bu t remains in use as a bench- 
mark (Ba stian et al.ll2010f) . the multisegmented power- 
law dKroupa et al.lll993l), and forms cl oser to log normal 
(|Miller fc Scaldll97ft ICharbrierll2003D . It is often quite 
difficult to disti nguish between multisegmen ted and log- 
normal variants (|Dabringhausen et al.l2008l) . Among the 
most direct measurements of the IMF, the evidence cur- 
rently points to a universal IMF that is a power law with 
index close to Salpeter's original value for stellar masses 
above a few solar m asses and log norm al or shallower 
for lower mass stars ()Bastian et al.l[2010f ), despite hints 
of variations in other, more extreme, environments or 
epochs that are observationally less accessible (for exam- 
ples s ee lRieke et al.lll993t !Davil2008t iNaravanan fe Pavel 
120121 ). In summary, the range of standard IMFs, let alone 
the potential for variations across different environments, 
results in stellar mass uncertainties of at least a factor 
of several, which is unsatisfactory when testing models 
whose purpose is to measure differences in galaxy stellar 
masses to comparable or better precision. 

To avoid the problems associated with testing these 
models against the observed properties of galaxies, with 
their complex stellar populations and ever present dark 
matter, we will test whether the models can reproduce 
the properties of local stellar clusters. This approach 
is not new. The challenge in implementing such a test 
is finding populations of clusters for which one can dy- 
namically measure the cluster masses and also sam- 
ple the relevant parameter ranges, principally cluster 
age, but also metallicity. To date, ~ 90% of the lo- 
cal clusters with kinematic measurements, from which 
masses can be derived, are old, l og(age [yrs]) > 10 
dMcLaughlin fc van der Marell[2005t iKruijssen fc Mieskd 
2009). The remainder tend to be extremely young and 
therefore subject to questions about whether t hey are dy- 
namically relaxed ()Goodwin fe Bastianl 12006). Because 
systematic effects can bias the estimated masses either 
high or low as a group, a population of clusters all at 
one age does not provide a strong test of the models. 
To test the evolutionary predictions of models, includ- 
ing the effects of dynamical evolution, investigators have 
turned to the pop ulations of clusters in galaxies outsid e 
the Local Group (|Reikuba et al.ll2007t lKruiissen|[2008l ). 
However, distance quickly diminishes our ability to mea- 
sure the internal properties of clusters, particularly age. 
Here, we seek, using Local Group clusters, to measure 
the evolution of the stellar mass-to-light ratio, Y*, for a 
set of clusters with well measured structural properties, 
ages, and chemical abundances, and to determine if mod- 
els of a simple stellar population with a universal IMF 
can produce a match to that evolution. 

A large compilat ion of Local Group stellar clus t er dat a 
was published by McLaug hlin fc van der Marel (|2005f) . 
The particular focus of that study was on the radial 
surface brightness profiles of the clusters and empha- 
sis was placed on homogenizing the data obtained from 



various sources, thereby making it useful for subsequent 
studies as well. In addition, when available, they col- 
lated velocity dispersions and produced dynamical mod- 
els. However, as we alluded to before, out of the 153 
clusters presented in that study, measurements of the in- 
ternal kinematics exist for only 57, and, of those, only 
6 are younger than 10 Gyrs old. They compared their 
dynamically measured T* to those derived from stellar 
population models, using bo th a variety of IMFs and 
two different evolution codes (|Fioc fc Rocca-Volmerangd 
119971 iBruzual fc Charlotll2003f ). Although the mean ra- 
tio of Y* from their dynamical models to that from their 
preferred stellar population models is 0.82 ± 0.07 (for- 
mally 2.5a discrepant with a ratio of 1), they do not 
stress this as a significant disagreement, presumably be- 
cause they appreciate that systemati c errors in their mass 
scale c ould be ~20%. Subsequently, Kruiisse n fc Mieskd 
(2009) alleviated the apparent discrepancy by showing 
that due to the internal dynamical evolution of the clus- 
ters, which causes the preferential loss of low mass stars, 
the measured values of Y* are expected to be lower than 
those calculated from stellar population models. Never- 
theless, even if the models, once one includes dynamical 
evolution, are in quantitative agreement for this set of 
clusters, this result provides support for the models only 
at a single age. 

To expand the age range of the sample of clusters 
with dynamically measured values of Y* and define a 
homogenous sample for study, we undertook to mea- 
sure the velocity dispersions of 22 Local Group stel- 
lar clusters. Here we present our measurements, us- 
ing integrated spectra, of the internal velocity disper- 
sions of 20 of those clusters that span 7 < log(age) < 
10.2 (for two we failed to obtain a measurement). Of 
those, we use the 18 that are also in the compilation of 
McLaug hlin fc van der Marell ([2005T) for our subsequent 
analysis. I n §2 we discuss the required data, b oth that 
drawn from lMc Laughlin fc van der M arel ( 2005|) and our 
own observations. We present the observational details, 
the data reduction, and our measurements of the internal 
velocity dispersions. In §3 we then convert these veloc- 
ity dispersion measurements to mass estimates. In §4, 
we present and discuss the behavior of Y* with age, and 
evaluate the possible effects of dynamical evolution and 
binary stars on these results. In §5 we focus on our pre- 
ferred hypothesis for the origin of the failure of the simple 
stellar population models to describe the behavior of Y* 
with age — the existence of at least two distinct IMFs. 
We summarize this study in §7. 

2. THE DATA 

Our spectroscopic data come from a set of observa- 
tions taken with the Las Campanas du Pont telescope 
(100-inch) and the Magellan Clay telescope (6.5m). All 
of the du Pont data and some of the Magellan data were 
obtained prior to the formulation of the current work and 
were meant primarily for an investig ation of the chemical 
abundance patterns within clusters (|Colucci et al.l[20TTl 
I2012T ) . For the bulk of the Magellan observations that are 
the core of this wor k, we selected clusters from the compi- 
lation provided by Mc Laughlin fc van der Marell (|2005h . 
The one exception is NGC 1718, which was observed 
(Nov. 2006) as part of the earlier work at Magellan. 
We use the compilation to select a range of clusters with 
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the necessary ancillary data (age, half light radius, lumi- 
nosity, modeled T* from stellar population models, and 
metallicity). From that list, we selected the clusters in 
Table Q] that were observed with Magellan on the ba- 
sis of their surface brightness within the half-light radius 
(which enables us to obtain high S/N, high dispersion 
spectra), with the additional requirement that we span 
the range of ages. Of the clusters that satisfy these crite- 
ria best, 6 already had measurements of the line-of-sight 
velocity dispersion, a, contained in the published compi- 
lation. Nevertheless, we observed those systems as well 
to obtain a homogenous, and homogeneously analyzed, 
set of velocity dispersion measurements. Finally, we also 
use the set of old clusters (age > 10 Gyr) with velocity 
dispersion measurements from the published compilation 
when we examine the question of dynamical evolution 
and its effect on Y» and discuss the possible origin of dis- 
tinct IMFs. When model data are used from the compila- 
tion we cho ose results ob tained using the Wilson mod- 
els (IWilsonl 119751 ). which IMcLaughlin fe van der Marell 
(2005) demonstrate are superior in fitting the radial sur- 
face brightness profiles of these clusters. Nevertheless, 
we tested whether our measurement of T* were affected 
by this choice and the results are discussed in §3. 

2.1. Ages, Photometry, and Stellar Population Models 

We adopt the age s, photometry, and stellar population 
results compiled by IMcLaughlin fc van der Marell (2005) 
without modification. The ages are derived from anal- 
yses of color-magnitude diagrams, compiled but not an- 
alyzed by those authors. The uncertainties are there- 
fore somewhat heterogeneous and, as always, dependent 
on the age of the cluster. We will broadly adopt a 
20% uncertainty in the age in our Figures. The pho- 
tometric measures, half-light radius, r/j, and mean sur- 
face brightness within the half-light radius, If,,, were cal- 
culate d and tabulated by IMcLaughlin fc van der Marell 
(2005) from their surface brightness fits and uncertain- 
ties are provided. The uncertainties in our dynami- 
cal estimates of T* are dominated by our uncertain- 
ties in the velocity dispersion er, partly because the 
fractional errors are larger on a than on any other of 
the required parameters and partly because the mass 
depen ds on a 2 . Finally, IMcLaughlin fc van der Marell 
(2005) present the results from a set of stellar popu- 
lation models. From the set they provide, we select 
those obtained with the iBruzual fc Chariot! (I2003f) m od- 
els and either a standard Salpeter I MF (ISa fpctcr 19551) or 
the often- used Chabrier disk IMF (|Charb ricr 2003J) and 
the lFiqc fc Rocca-Volmeran gc (1997) (PEGASE) models 
and a iKroupa et all ()1993D IMF. We adopt the uncer- 
tainties quoted in the compilation, although stochastic 
effects (how well sparsely populated features in the CMD 
are sampled) can lead to significant deviations from the 
expected photometric properties. For our clusters, which 
have 3.5 x 10 4 < M < 5.0 x 10 5 A/ Q , as we will show be- 
low, the role of stochasticity in colors is potentially signif- 
icant, although it is drama tically larger for M < 10 4 M Q 
([Pci escu fc Hansonll2010al ). and" some scatter in the total 
luminosity is also expected. 

We will also go beyond the model results provided 
by IMcLaughlin fc van der Marell (120051) by utilizing the 
models provided by I Anders et al.l (|2009l ). which address 
the issue of the dynamical evolution of stellar clusters. 



We discuss these in detail in ij4.ll 

2.2. Integrated Spectra 

We obtained our spe ctra using the MIKE spectrograph 
(jBernstein et al.l I2003D at the Magellan II (Clay) tele- 
scope during Nov. 2006, Aug. 2009, Feb. 2011, and Oct. 
2011, and with the echelle spectrograph at the 100-inch 
at Las Campanas (du Pont) during Dec. 2000 and Jan. 
2001. 

The MIKE spectrograph provides blue and red chan- 
nel spectra over the entire optical window. We use the 
data from the red side, in which the most prominent ab- 
sorption lines with high S/N lie, and confine ourselves 
to shortward of 7000 A to avoid strong atmospheric fea- 
tures. We used the 0.75 arcsec slit to avoid degrading 
resolution further, which is already near the limit for 
clusters that can possibly have a ~ 2 km s _1 . The expo- 
sure times vary depending on the surface brightness and 
are given in Table [T] 

To obtain integrated spectra of the clusters, we uti- 
lize the same spectr oscopic drift technique described by 
iColucci et al.l (|2011[) . To summarize, we set the telescope 
in motion to raster scan the slit across the cluster dur- 
ing the exposure, defining both the angular length and 
height of the raster (both set to the same number). The 
exposure time then sets the rate of the scan, such that 
the full scan is completed within the allotted exposure 
time. We chose from only two different values of the 
scan sizes (10" x 10"or 30" x 30"), aiming for the closest 
match to the half light radius of the cluster. Because 
MIKE does not have an instrument rotator, these scans 
are then further complicated by field rotation, which is 
unaccounted for. Multiple exposures further help homog- 
enize the sampling of the central region. Our goal is to 
obtain a spectrum that is representative of the central 
region of each cluster rather than one of a tightly speci- 
fied region. The slit is only 5 arcsec long, so these raster 
scans have no clear sky. We do not attempt sky sub- 
traction. Solar spectral features, from moonlight, are at 
significantly different velocities and are also quite weak 
as the two key runs (the last two) were both in dark time 
(moon illumination < 30%). 

We reduced the MIKE data with the MIKE DR 
pipeline ( R. A. Bernstein, in p reparation, available at 
http://www.ucolick.org/ ^xavier/Profession.html). Be- 
cause of the scanning technique, the cluster light fills the 
slit, so that a modification to the pipeline was developed 
that allows boxcar rather than optimal source extraction. 
The data obtained in Aug., Feb. and Oct. 2011 gener- 
ally have shorter exposure times, fewer frames and more 
intrinsic emission lines from the youngest clusters, which 
render the standard cosmic ray rejections in the MIKE 
DR and IRAFQ packages ineffective for producing the 
combined spectra. For these data, we remove the cosmic 
ray and emission spikes in the individual frames using 
a 10-sigma threshold above a continuum fit and median 
combine the cleaned frames using the IRAF scombine 
routine. Again, we remove the blaze function using a 
trace of the continuum flux of a G-type star. 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 
Stellar Cluster Data 
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A second set of data come from observations with the 
du Pont telescope and its echelle spectrograph. Those 
spectra have a wavelength coverage of approximately 
3700 to 7800 A, with declining sensitivity and spectral 
resolution toward the blue end. These spectra were also 
obtained using the scanning technique. At the du Pont 
telescope scanning was implemented using a modification 
of the telescope guider program provided by S. Shectman 
(see Mc William & Bernstein, 2008). The echelle slit is 
l"x 4", allowing uniform coverage of a 12 x 12 arcsec 2 
or 8 x 8 arcsec 2 high-surface brightness region of the 
cluster. We took multiple exposures to homogenize the 
scanned cluster region and for cosmic ray removal. We 
reduced the spectra using standard IRAF routines (see 
Colucci et al. 2011), including the scattered light sub- 
traction described in detail in McWilliam & Bernstein 
(2008). We combined the extracted spectra using the 
IRAF scombine routine with the crreject algorithm to 
eliminate cosmic ray events. Finally, we remove the blaze 
function using a trace of the continuum flux of a G-type 
star. 

2.3. Measuring Velocity Dispersions 

Various techniques have been developed and applied 
to the problem of measuring velo city dispersion s from 
integrated spectra (for example , iSargent et al.l 119771; 
Benderl [l990t iRix fc Whltel Il99l iKuijken fc Merrifieldl 
19931 : IWinsall &: Freeman! 11993) . The appropriate ap- 
proach depends on the particulars of the scientific situa- 
tion and data quality. In our case, we have good S/N for 
most of our clusters, and therefore the luxury of not hav- 
ing to rely on a method that utilizes the entire spectrum 
at one time so as to increase the S/N, but which incurs 



the cost of obscuring potential systematic errors. Specif- 
ically, we have enough S/N that one trustworthy line 
would be sufficient to provide a results with the desired 
precision. As such, our primary concern is systematic 
errors. 

Systematic errors in this type of measurement can arise 
from a variety of sources including differential instrumen- 
tal broadening across the spectra, template mismatch, 
errors in wavelength calibration, and errors in the defini- 
tion of the continuum level. These are all much easier to 
identify if one can work with individual absorption lines 
rather than with a vaguely weighted mean from the com- 
bination of all lines. Some of these errors can be amelio- 
rated by fitting on a case-by-case basis, for example one 
can let the radial velocity float in the fits of individual 
lines to mitigate against wavelength calibration errors, 
minimizing the risk that an inflated a will be inferred as 
the algorithm attempts to better fit a number of lines 
with a single radial velocity. Other errors, such as those 
that arise for lines near the ends of orders, will result 
in "problem lines" being identified as outliers in the set 
of measurements. For these reasons, we choose to fit as 
many of the absorption lines as we can individually and 
then use the entire set of measurements to identify the 
most likely value of the velocity dispersion. 

Our selection of lines to be fit is based on a visual in- 
spection of a long list of lines, originally identified in our 
template star spectra. We reject any lines in the object 
spectra that are either clearly blends or suffer some other 
complication and those that do not dip below 0.85 of the 
local continuum level. After fitting, we reject lines that 
do not produce an acceptable fit, where acceptable is de- 
fined by x 2 < 2.3. To calculate % 2 we adopt a per pixel 
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uncertainty determined from the fluctuations about a fiat 
continuum in line-free areas of the spectrum. However, 
we adopt the same uncertainty value for the full spectral 
range for any given cluster. Our results are not highly 
sensitive to this value because we only use these % 2 val- 
ues to remove questionable fits from further considera- 
tion, and visual inspection confirms that those lines that 
have been rejected by this criteria are poorly fit. 



NGC2002:HR0805 




Wavelength 



Fig. 1. — Example absorption line fits drawn from data for NGC 
2002, using stellar template HR 805. Each line is fit independently 
using a Gaussian convolution kernel and scaling for both slightly 
different radial velocities and equivalent width. Best fits are shown 
with the smooth line (red). These are a subset of the absorption 
lines available for NGC 2002. 

Some of the subtleties involved in selecting among the 
various methods in the literature to measure a lie in the 
possible nature of the underlying line-of-sight velocity 
distribution (LOSVD), which can be quite complex in 
galaxies. It is our expectation that in clusters, with 
relatively short dynamical times, the LOSVD is nearly 
Gaussian. As such, our fitting involves the convolution 
of a Gaussian, which we adopt as the description of the 
LOSVD, of specified width with the rest frame template 
spectrum. For our MIKE cluster spectra, we use two 
template stars and later compare results (we find no sig- 
nificant differences in the inferred a). The convolved 
spectrum is then redshifted by a specified value corre- 
sponding to a selected radial velocity, v r , and binned 
into the same wavelength bins as the object spectrum 
using cubic splines. Finally, we also explore a range of 
equivalent width renormalizations so that we can match 
the strength of the line in the template to that in the ob- 
ject spectrum. The renormalization is done by subtract- 
ing 1, the value of the normalized continuum, from the 
spectrum, multiplying the spectrum by a specified fac- 
tor, and adding the 1 back to the spectrum. We sample 
the multidimensional parameter space (cr, v r , normaliza- 
tion) uniformly and in an unbiased way, calculate x 2 for 
±3 pixels about the line center (in a few cases where the 
spectral lines are narrow we set this to ±2 pixels) because 
we do not want the selection of the best fit parameters 
influenced by nearby spectral lines or slight continuum 
mismatches, and identify the best fit parameters. The fo- 



cus on the core of the line minimized various systematic 
technical problems, such as fluctuations in the contin- 
uum, but also addresses some external sources of uncer- 
tainty, such as the role of binaries (§4) and interlopers, 
by placing less weight on the wings of the line-of-sight ve- 
locity distribution. Examples of fits for absorption lines 
in the spectra of NGC 2002 are provided in Figured] 

We correct for instrumental broadening using the tem- 
plate star spectra. Each of the absorption lines we iden- 
tify in the template spectrum is fitted with a Gaussian. 
Assuming that the intrinsic line widths are well below 
the instrumental broadening, which is in the range of 2 
km s _1 , then the fitted Gaussian widths as a function 
of wavelength describe the instrumental broadening. We 
do indeed find a strong correlation between the Gaussian 
width in km s _1 and wavelength, as expected from in- 
strumental broadening. We fit a low order polynomial 
to that relation and use it to estimate the instrumental 
broadening at any wavelength. This procedure is done 
independently for each run and, of course, each of the 
two instruments used. The inferred instrumental broad- 
ening is added in quadrature to the best fit value of a 
derived above because we have fit a template line that 
has already been broadened. 

Once the line fitting is complete, the question becomes 
how to use the values derived for the set of individual 
absorption lines to best derive the value of a. The use 
of individual lines, even when visually vetted, does not 
eliminate the potential for systematic effects that inflate 
the widths (such as blends and focus errors). The in- 
ternal errors provided by the fits themselves are usually 
optimistic given the high S/N of these data. We eval- 
uate the final value of a for each cluster by calculating 
the average of all the measurements and % 2 /N , where N 
is the number of data points. We set the uncertainty in 
any individual measurement by requiring x 2 /A = 1. To 
downweight measurements that are inflated by blends or 
focus errors, we set upward uncertainties to be X times 
larger than downward ones. We will discuss the selec- 
tion of X shortly. The uncertainty on our final "mean" 
a is derived by identifying the range of cr's that generate 
X 2 — Xmin < 2.71, corresponding to the 90% confidence 
interval (see Table [TJ. 

Different weighting factors X will result in different 
values of a. As such, without further constraints on 
X, this parameter is itself a source of systematic un- 
certainty, and simply reflects the uncertainty associated 
with unknown blends, focus errors, etc. To the degree 
that changing X results in similar proportional shifts to 
all the cr's, this effect only results in global shifts in the 
values of Y», whose normalization is uncertain for other 
reasons as well. Nevertheless, for guidance we exam- 
ine the internal uncertainty on the individual fitted lines 
produced by our fitting algorithm, which have upward 
uncertainties that are between 1.5 and 3.5 times larger 
than downward uncertainties. For our calculations, we 
therefore select a value within that range, X = 3, and 
note that choosing X anywhere in this range produces 
resulting cr's that lie within the quoted uncertainties in 
Table [TJ 

In certain panels in Figure [2J pairs of points are con- 
nected with lines. These points and lines represent two 
measurements of the same absorption line using different 
template stars. We find that on average, over the twelve 
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Fig. 2. — Results from the velocity dispersion determination. In each panel we present the results from our measurement of velocity 
dispersions using individual absorption lines that passed our selection criteria for a different cluster. In cases where we also compare results 
using different template stars, we use filled and open circles and connect the results with vertical lines. We find no significant systematic 
difference in <r arising from our choice of template. The best fit dispersion value for each cluster is shown as the solid line, with the grey area 
surrounding it demarcating our 90% confidence interval on this mean value. The dashed lines show the uncertainty, derived as explained 
in the text, for an individual measurement. 



clusters for which we used two templates, the velocities 
of individual lines are biased by 0.12 km s _1 from one 
template star to another — a value smaller than any of 
the internal uncertainties we quote. Therefore, we con- 
clude that template selection does not have a dominant 
effect on our uncertainties. 

For the subset of our clusters that have been ob- 
served previously, we compare our a measurements to 
pre-existing values. First, one of us (JC) has indepen- 
dently (although using the same spectra) obtained es- 
timates of a using a cross-correlation analysis. This 
provides an internal test of our methodology. Sec- 
ond, fo r an external test, we turn to th e compila- 
tion of McLaughli n fc van der Marell (|2005l ) , which in- 
cludes literature values of a for six of our clusters, and 
other sources tha t include two more (|Mevlan et alJll99ll : 
lLane et "all I2010I ). Both internal and external tests are 



presented in Figure [3] The external comparison is pos- 
itive in that 5 of 8 (63%) of the measurements agree to 
la and 7 of 8 (88%) ag ree to 2a. For the one significant 
outlier NGC 1866, the lLane et al.1 ((2010) study can be 
used to demonstrate the sensitivity of that particular re- 
sult, obtained using individually measured stars, to the 
inclusion or exclusion of single star (dashed line in Figure 

El). 

There are two clusters, NGC 458 and 1860, for which 
we were unable to identify any suitable lines for mea- 
surements. This is somewhat surprising given exposure 
times that are in the range of those obtained for other 
clusters. NGC 1860 is particularly large, nearly twice as 
large as the next nearest cluster in size, and has a corre- 
spondingly low mean surface brightness, which may have 
just been too low for these observations. The situation 
for NGC 458 is a bit more puzzling, but we may have 
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simply been unfortunate in the region we scanned. 
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Fig. 3. — A comparison of our measurement of a to two other 
measurements. First, we compare to determinations using our data 
but a different technique (cross-correlation) in the filled circles. 
Second, we compare to existing m easurements in the literature 
HMcLaughlin fc van "der Marcl 2005) in the filled triangles. The 
line represents the 1:1 relation. 



3. DETERMINING MASSES 

The single-epoch dynamical determination of the mass 
of an astronomical system always references the Virial 
Theorem. In practice, the different methods that have 
been developed are all attempts to deal with the un- 
known numerical factors arising from the integration of 
the mass distribution and kinematics required when cal- 
culating the terms in the Virial Theorem. In certain 
cases, where the number of independent tracer parti- 
cles with measured velocities is small, there is a con- 
cern for numerical stabili ty of the mass estimator (cf. 
IBahcall fe Tremaind [l9811i). but usually the concern is 
how to properly weight the kinematic term, related to 
a - , given the unknown distribution of tracer particle or- 
bits. Theoretically, the difference between the mass in- 
ferred for a system that has a tangential distribution 
of orbits rather than a radial one can be as large as a 
factor of three if the incorrect orbital anisotropy is as- 
sumed. There are observational ways to constrain the 
orbital anisotropy using higher or der measures of the 
line of sight velocity distribution (seel Winsall fe Freeman] 
119931: iKuiiken fc Merrifield 1993), but such analyses are 
beyond the scope of the current work. 

A surprising work -around to this pro ble m has been 
identified recently bv lWalker etaLl (|2009[ ) and lWolf et all 
(|2010t) . They find that the mass enclosed at the half 
light radius of stellar systems, and only the half light 
radius, is nearly insensitive to such details as the or- 
bital anisotropy or the radial profile of the tracer pop- 
ulation. They have gone on to use this finding to es- 
timate the masses of dwa rf spheroidal galaxies in par- 
ticular (Walk er et al.ll2009f). but also of other spheroidal 
stellar systems (|Wolf et^!1[20T0l). By examining a range 
of dynamical models I Walker et al.l (|2009f) found a robust 



estimator for Mh to be Mh — 580r/i<7 2 , where M is in 
solar masses, rh is in parsecs, and a is in km s -1 . 

Such an estimator is even more trustworthy if verified 
empirically. Of particular importance is whether the es- 
timator is valid over the full range of system masses and 
the degree to which mass estimates for individual sys- 
tems fluctuate about the me an value descr i bed b y the 
estimator. Interestingly, the Walk er et al.l (|2009l ) esti- 
mator, w hen rewritten in the style of the Fundamental 
Manifold ( Zaritsky et al. 2008), is equivalent to 

log r h = 2 log a - log I h - log Y/j - 0.73, (1) 

where now is given in kpc, a is still in km s _1 , Ih, 
which is the mean surface brightness within rh is in solar 
luminosities per sq. parsec, and is the mass-to-light 
ratio within r h- By measuring rh, o~, and , Ih, one can 
solve for T h . IZaritskv et all (|2008l 120111) identified an 
almost ide ntical empirical scali ng relationship, indepen- 
dent of the lWalker et al.l (|2009D study, that works for all 
stellar systems, ranging from the most massive galaxies 
to globular clusters. They expressed their relationship as 

log r h = 2 log a - log I h - log T h - 0.75. (2) 

where the zero point, the 0.75 term, comes from placing 
galaxies with indepen dent mass estimates fr om detailed 
dynamical modeling (|Cappellari et al.ll2006l ) on the re- 
lationship. As is already evident, this empirical calibra- 
tion produces an almost exact match in the normaliza- 
tion (0.73 v s 0.75) to that prov ided by the theoretical 
modeling of I Walker et al.l ((2009). The empirical results 
verify that there is little scatter (~ 0.1 dex) about this 
relationship for objects ranging from globular clusters to 
massive elliptical galaxies. Using Equation [51 we evalu- 
ate Th for our set of 20 clusters and present the results in 
Table [1] Masses can be calculated using Th and the to- 
tal luminosities, which we also provide in the Table. For 
systems without dark matter, which we presume includes 
these clusters, Y* = Th- All photometric quantities are 
presented for the V band. 

The uncertainties in Y* are calculated using only the 
uncertainty in a, as the internal uncertainties on the 
other parameters are proportionally much smaller. How- 
ever, we test for the possible systematic effect in 
and Ih due to our choice of the Wilson model fits by 
comparing estimates of Y* using t he values given by 
McLaughlin & van der Marcl (2005) for their alternate 
King and power-law fits. We find that with the excep- 
tion of two clusters (NGC1856 and NGC2100), where the 
power-law fit gives significantly different estimates of rh 
and Ih, all estimates of Y^ are within 2a and 80% are 
within lcr. For the two discrepant clusters, the power- 
law fit to NGC1856 results in an unphysically large value 
of rh, and hence Y„, and in both cases the results of the 
King and Wilson models agree within the uncertainties. 

3.1. Comparison to Other Estimates ofTh 

We compare our estimates of Th using Equat ion[2]with 
those obtained by iMcLaughlin fc van der Marell |2005) . 
Using the published values of a, but our mass estimator, 
we obtain va lues of Y* that tend to be low e r than those 
obtained by IMcLaughlin fc van der Marell (|2005l ) by < 
30% (Figure [4]). However, global shifts of this magni- 
tude are expected when using different estimators due to 
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the difficulty in determining the zero point calibration 
of such estimators. More importantly, the clustering of 
points about a fixed ratio of the two estimators shows 
that the correlation between the two is excellent. If we 
were to recalibrate Equation[2]to give consistent answers, 
we would change the calibration constant in the Equation 
from 0.75 to 0.60. We will leave the consta nt as is, partly- 
due to the agreement with the results from I Walker et al] 
( 2009), but appreciate that there is a degree of freedom 
available in the overall normalization of the values of TV 
That freedom is a global shift in T/, rather than one that 
depends on itself. 
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F ig. 4. — The dependency of the ratio of estimates 
from IMcLaughlin fc van der Marell (f2005T) (THEIRS) to our own 
(OURS) on central projected mass density. The bottom portion of 
the Figure illustrates how th e our estimate of Y^ is well b e haved 
relative to that provided by McLaughlin & van dcr Marel (2005) 
for central projected mass densities that are < 10 4 Mq. There is 
a systematic offset where our values tend to be lower by ~ 30% 
(the dashed lines shows a ratio of 0.7 for reference). Above that 
value there is large scatter, with an implication that our method 
could produce significant overestimates of Tj, for some large val- 
ues of the central mass surface density. The top portion of the 
Figure illustrates the distribution of our sample with respect to 
central projected mass density. Only three of our clusters lie in 
the suspect region, although those three are all old (log(agc) > 10) 
clusters and our estimates of T* for these are low. 

Examining the differences in detail (Figure [4]), we find 
that there is a marked connection between the difference 
in the values of for individual clusters and their cen- 
tral projected mass density (and related quantities like 
volume density and phase space density). Below central 
mass surface densities of 10 4 M Q pc~ 2 , the two mass es- 
timators track each other well, modulo this < 30% offset. 
For central mass densities greater than this, the scatter 
increases significantly. This does not, in itself, suggest 
which of the two approaches should be preferred, but sug- 
gests some caution in treating estimates of for high 
density clusters. In our sample, only three clusters (NGC 
1916, 2005, and 2019; see red, open symbols in the plot) 
lie above this threshold. We find no connection between 
the differenc es in mass estimates and other qua ntities 
provided by IMcLaughlin fc van der Marell (|2005| ). such 
as relaxation time at the half mass radius or total mass. 



The d iscrepancy, if the IMcLaughlin fe van der Marell 
(2005) estimates are the ones to be trusted, is in the 
sense that we would be overestimating for some frac- 
tion of high central density clusters. The three clusters in 
our sample that are potentially at risk are all old clusters 
(log(age [yrs1)>10), but our estimate s of and those of 
IMcLaughlin fc van der Marell (|2005[) agree to within 5% 
(after adjusting for the 30% mean offset). 

4. STELLAR MASS-TO-LIGHT RATIOS 

Using the derived values of Y/j and assuming that the 
clusters are devoid of exotic dark matter, we show our 
principal empirical result, the relationship between T* 
and age, in Figure [S] Here we see a definite trend, where 
Y* rises as we consider older and older clusters up to 
ages of a few Gyr, and then drops significantly for the 
oldest clusters. While this drop runs counter to the naive 
expectation that clusters should continue to fade with 
age and therefore that Y* should continue to rise, the 
evaporation of low mass stars due to two-body relaxation 
will alter this expectation, particularly at the oldest ages, 
as discussed in §4.1. 








log(age) [yrs] 



Fig. 5. — Stellar mass-to-light ratio, Y* versus age. We plot 
the dynamical estimates of Y* for our homogeneous sample of 18 
stellar clusters. 



Before continuing to explore all of this further, we note 
that the youngest two clusters in ou r sample are close to 
or bel ow the 50 Myr age for which iGoodwin fc Bastianl 
(2006) argue that dynamical mass estimates are unreli- 
able because clusters are still undergoing violent relax- 
ation. They show that dynamical masses can be inflated 
by factors of several, although these effects go away for 
log(age [yrs])> 8. We therefore place no weight on dis- 
agreements found between observations and models at 
the youngest ages (log(age [yrs]) < 8. 

We compare our results to models of simple stel- 
lar populations (Figure El) . Specifically, we begin with 
the re sults presented by IMcLaughlin fc van der Marell 
(2005) using either thelBruzual fc Chariot! (I2003D or PE- 
GA SE dFioc fc Rocca-Volmerajigdll997fr ' models and ei- 
ther iSalpeterl (|1955l ). iCharbrierl (|2003l) . or iKroupa et al.l 
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Fig. 6. — Stellar mass-to-light ratio, Y*, versus age compared to 
model prediction. To Figure [5] we have now added the calculated 
val ues of T« for different st ellar population models as evaluated 
by IMcL aughlin & van dcr Marcl (2005]) . The upper (green) curve 
represents the re sults obtained u sing the Bruzual & Chariot; (120031 ') 
algorithm and a Salpctcr (1955) IMF. The i ntermediate (red) rep- 
resents results using the same algorithm by a Charbricr (2003) disk 
IMF. Finally, the lowest cu r ve rep resents results obta ined using the 
Fioc & Rocca-Volmcrangc (1997) algorithm and the Kroupa et al. 
(119931 ) IMF. All three fail in quantitative detail to fit the data. 
The clusters with log(age) < 8 may not be relaxed and are not 
given any weight in this comparison. Dynamical effects that may 
affect the older clusters are discussed in the text, but found to be 
insufficient to reconcile the data and models. 



(1993) IMFs. Estimates of T* were calculated using the 
metallicities listed in Table [TJ and so in theory account 
for the different metallicities of the old and young clus- 
ters. The models clearly straddle the data (FigureE]) but 
fail both by underestimating T* for the younger clusters 
and overestimating it for the older clusters. Simple mul- 
tiplicative rescaling of T* is allowed because our normal- 
ization (Equation^ could be questioned — for example, 
our calibration based on elliptical galaxies might not be 
an exact match for stellar clusters — but a simple multi- 
plicative shift of the T» measurements cannot address an 
underprediction at young ages and an ovcrprediction at 
older ages. The only way to reconcile the data and mod- 
els is to posit that either the observations or the stellar 
population models have a systematic error that depends 
on age (or a related parameters such as metallicity). 

One might suspect the results because there are only 
a handful of points at the oldest ages. Perhaps these are 
an omalous among older clusters . The sample provided 
by IMcLaughlin fe van d cr Marcl (|2005D contains 51 old 
clusters (log(age [yrs]) > 10) and the median of the T* 
distribution for those clusters, when estimating Y* with 
Equation [51 is 1.2, or about 40% larger than the aver- 
age T* for our five old systems (0.85). This difference 
includes the 30% offset discussed in the previous section 
that is the result of a zero point mass offset, so the addi- 
tional factor is small. We conclude that while our clus- 
ters may have somewhat unusually low values of T„, the 
likely mean value still lies well below the peak T* seen 
at intermediate ages, and certainly does not lie on the 
rising extrapolation of the relation seen for the younger 



clusters. 

There are, of course, physical reasons why the results 
may be distorted. First, the older clusters are all of sig- 
nificant lower chemical abundance than the younger clus- 
ters (see Tabled]). This is certainly a source of concern, 
although we are comparing the values of Y* to those pre- 
dicted by models that account for metallicity differences. 
Second, these clusters could, for some unknown reason, 
not satisfy the mass estimator of Equation [2j Aside from 
the high central projected surface brightnesses of three 
of these, there is no striking difference in structural pa- 
rameters. They tend to be among the physically most 
compact, most concentrated ones, but not exclusively so 
(NGC 121 and NGC 4590, fall within the range of the 
majority of the clusters). Furthermore, several of these 
(NGC 1916, 2005, 2019, and 4590) are also among the 
clusters studied bv IMcLaughlin fe van der Marell (2005) 
and the mean ratio of our T* estimates to theirs is 0.6. 
Recall that for the entire sample of old clusters this ra- 
tio is ~0.7, so we see no unusual behavior in the esti- 
mated T» values for our set of old clusters beyond the 
increased scatter at high surface brightnesses that was 
discussed previously. If anyth i ng, th e comparison to the 
IMcLaughlin fe van der M arcl (2005) results, which high- 
lighted the issue with high surface brightness clusters, 
suggests that we are likely to be overestimating T* for 
this set of clusters. Disregarding this comparison and 
hypothesizing a factor of two upward scaling error in our 
Y* values would allow our cluster T* estimates to come 
to better agreement with the evolutionary models that 
use Chabrier or Kroupa IMFs, but would exacerbate the 
underprediction of Y* for the intermediate age clusters. 
We find no evident reason why the data for the older clus- 
ters in our sample should be singled out as susceptible 
to calibration biases in Y* by factors > 2. 

Alternatively, one could be suspicious of the intermedi- 
ate age clusters. However, the values of <r, rh, and Iu for 
these clusters are roughly within the value ranges for the 
other clusters (Figure[7|), with the possible exception that 
their surface brightnesses are somewhat lower. The latter 
is due to the lack of intermediate age clusters available for 
study and our selection for the highest surface brightness 
clusters available. Even so, this range of surface bright- 
nesses is not outsid e of the range of our comparison to the 
modeling results of IMcLaughlin fc van der Marell (2005) 
(Figure 0]) and so we do not see this as a source of con- 
cern. 

The only way to reconcile the data and models is to 
posit that either the observations or the stellar popu- 
lation models have a systematic error that depends on 
age (or a related parameters such as metallicity). We 
have just argued against systematic errors in our sam- 
ple selection or measurements, and so we are left with 
the possibility of age-dependent problems in the applica- 
tion of simple stellar population or dynamical models to 
predict Y*. We address two potential reasons for such 
failures below. 

4.1. Dynamical Evolution of Clusters and Its Impact on 

Y 

When using stellar clusters, particularly as a test of 
stellar population models, one needs to account for dy- 
namical effects that lead to the preferential loss of low- 
mass stars and alter the simple predictions for Y* shown 
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Fig. 7. — Cluster structural properties as a function of age. The 
old and intermediate age clusters do not have strikingly different 
structural properties. 



in Figure |U Clusters can be significantly affected by 
mass loss driven by gas removal at early times (for ex- 
amples see lHilldfl^80TlGoodwin fc Bastianll2006f ). and by 
stellar evolution in a tidal field, gravitational perturba- 
tions due to passing molecular clouds, ejections during 
binary interactions, and e vaporation o f low mass stars 
via two-body interactions (|Spitzerlll958t lKruiissenll2008fl 
on longer timescales. Modern treatments of the evolution 
of an isolated cl u ster or one in a tidal field are presented 
by iGieles et ail (|2010al 1201 If ). We concentrate on the 
evaporation of low-mass stars via two-body interactions, 
which has the strongest impact on the oldest clusters, 
to determine the degree to which the measured values 
of have been lowered relative to the simple stellar 
populations expectations. 

Low-mass clusters are disrupted earlier than high-mass 
clusters of similar density due to internal two-body re- 
laxation. The evaporation rates depend on a number 
of poorly constrained parameters, such as the internal 
density profile and tidal field strength. It is therefore 
difficult to predict the evaporation rate for any partic- 
ular cluster. Instead, the luminosity and mass function 
of globular clusters provides constraints on the evapora- 
tion rate for populations of clusters. The mass function 
of globular clusters in the Milky Way has a peak near 
« 1 - 2 x 10 5 M Q . Fall k Zhang (2001) dynamically 
evolved simulated globular cluster systems and showed 
that this peak results from the earlier disruption of low- 
mass clusters relative to high-mass clusters of similar 
density, where the mass of each cluster is depleted ap- 
proximately linearly with time, with an evaporation rate 
fi ev w 1 — 2xl0~ 5 Mq yr^ 1 . The mass function of globu- 
lar clusters in the LMC has a similar shape to those in the 
Milky Way, although the uncertainties on this shape are 
larger because of lower numbers (Chandar et al. 2010). 
Nevertheless, this result suggests that the ancient glob- 
ulars in our LMC sample have experienced similar evap- 
oration rates, on average, as those in the Milky Way. 

Various authors have studied the impact of relaxation 
driven evaporation on T* in star clusters. For example, 
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Fig. 8. — Correcting for the effects of dynamical evolution on X*. 
In the upper p anel w e show the results from applying the models of 
lAnders ct al. (2009), which address the issue of dynamical evolu- 
tion of the clusters. The correction is in the sense of accounting for 
lost low mass stars from clusters so as to provide a fair comparison 
to simple stellar population models. We plot the ratio of T* ob- 
tained from a model with effectively no evolution, setting fcj = 100 
Gyr, to ones where we have taken the lowest, self-consistent val- 
ues of ti (see text for details). The dynamically derived values of 
T* would need to be multiplied by the plotted value to obtain an 
estimate of Y„ for a population with no lost stars. In the lower 
panel, we apply that correction and find that although there is 
an increase in the values of T* for the older clusters, the quali- 
tative nature of the behavior of Y, with age remains unchanged 
and will still not match the expectations drawn from simple stellar 
population models. 



Lamers et al. (2005a) and Anders et al. (2009) used N- 
body simulations of star clusters to predict total cluster 
lifetimes and the evolution of Y* for clusters with dif- 
ferent stellar mass- loss rates. Relative to standard mod- 
els, low-mass stars, which have a higher Y than that 
of the typical cluster star, are lost, but the fraction of 
non-luminous stellar remnants increases and these two 
effects partially cancel. Anders et al. (2009) present cal- 
culations for the disruption time for clusters, td, which 
depends on the cluster mass and local (external) density, 
with a mass dependence of td oc M 62 . This relation- 
ship is typically normalized by specifying the time at 
which a 10 4 M cluster loses 95% of its mass, ti. The 
mass dependence can be converted into a mass loss rate 
dM/dt oc M 38 , which is shallower than the relation 
given by classical evaporation and gives a poor fit to the 
overall shape of the mass fu nction of globular cl usters in 
M87 and in the Milky Way (fWaters et al.ll2006f) . partic- 
ularly at the low mass end. Nevertheless, these models 
can reproduce the peak of the mass function. In prin- 
ciple, each cluster has its own corresponding value of t\ 
depending on its internal density, orbit, and tidal field 
strength, although we have insufficient information to 
evaluate t± on an cluster-by-cluster basis, and a single 
value of ti is usually assumed for a set of clusters. 

To evaluate the ma gnitude of th e effec t on T„, we 
choose from among the I Anders e t al. (2009) models for a 
selected value of £4 (see below) and the nearest values of 
(Fe/H) and age for the particular cluster. In Figure [5] we 
show the ratio between the predicted values of Y* from 
these models to those obtained from models with effec- 
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tively no evolution (£4 ^100 Gyr). To select values of 
£4 we adopted values from the liter ature. In the SMC, U 
is estimated to be roug hly 10 Gyr (ILamers et al.ll2005bD . 
while in the loc al neighborhood it is estimated to be be- 
tween 1.3 Gvr ([Lamers et al.ll2005af) and ~5 to 6 Gyr 
(|Baumgardt fc Makinol 120031 : lLamers et al.l l2005b| ) . We 
adopt £4 = 10 Gyr for clusters in the SMC and then scale 
by mass for our SMC clusters. We adopt t 4 = 5 Gyr for 
the remaining clusters, which is the lowest possible value 
of ti that produces self-consistent results in the sense that 
none of our clusters are older than the inferred dissolu- 
tion time. The results shown in Figure [H] demonstrate 
that the effect is at most a factor of two increase in T* 
although usually significantly less. Hereafter, our values 
of T* are the corrected values. A comparison between 
the present-day stellar mass functions measured for 27 
old globular clusters and simulations of cluster evolution 
suggest that the present-day mass functions are consis- 
tent with expectatio ns from two-body relaxation for a 
given universal IMF (jLeigh et al.ll2012h . 

There is, of course, uncertainty in the selection of £4. 
However, there are two reasons why we conclude that 
dynamical evolution cannot reconcile the models to the 
data even if one is allowed to change £4 on a cluster-by- 
cluster basis. First, the effects of dynamical evolution 
are at best only modest (~ 2) and reconciling the data 
to the better fitting Salpeter model (Figure [5]) requires 
significantly larger corrections for the entire sample of old 
clusters. Second, despite the choice of £4 within reason, 
some of these old clusters are sufficiently massive that 
they are impervious to the effects of dynamical evolution. 

5. BINARIES 

If the old clusters are behaving as expected, then the 
problem might lie with the intermediate age clusters. If 
so, one possibility is that their velocity dispersions are in- 
flated by binaries. Contamination by binaries, which can 
have orbital velocities > er, have been a long running con- 
cern in the measurements of a for dwarf spheroidal galax- 
ies because those mea surements imply tremendous quan- 
tities of dark matter (lOlszcwski et al. | ll996t [Minor et all 
120101 ) and for young stellar clusters (|Gieles et alJl2010b| ). 
On the other hand, the issue has not raised much con- 
cern in the analysis of globular clusters because the in- 
ferred masses for the old clusters are in moderate agree- 
ment with expectations drawn from stellar models. Di- 
rect investig ation of the binary fraction in a well-studied 
old cluster (|Prvor et al.l 119881 ) has confirmed a low bi- 
nary fraction, b ut the fractions are found to be larger in 
young clusters dPortegies Zwart et al.l l2010h and in low 
mass clusters (|Milone et al.ll2012f) . We need an estimate 
of the magnitude of the effect binaries might have on the 
measured cr's of our clusters. 

iMinor et~all ()2010l ) did extensive modeling to address 
this issue for low mass dwarf galaxies. Except for possi- 
ble internal evolution in clusters, such as the destruction 
of binaries and mass segregation, the results of these cal- 
culations should be applicable here. They find that for 
a system with an intrinsic value of a = 4 km s , which 
is appropriate for our intermediate age clusters, bina- 
ries are likely to bias a upward by between 10 and 20% 
depending on the binary fraction. Again, this model- 
ing depends on the various properties of binary stars in 
addition to the fraction, but the fraction may actually 



be lower in clusters than adopted as typical elsewhere 
due to the increased dis ruption of binaries in clusters 
(|Spitzer fc Math eiu 1983). Even though this level of bias 
is comparable to or slightly larger than the uncertain- 
ties we quote, it is small compared to the size of the 
effect that concerns us, which is at least a factor of ~ 2 
larger. The effect should be further reduced in our mea- 
surements relative to their calculation because we place 
little weight on the wings of the LOSVD by fitting the 
line cores (this also helps mitigate against interloping 
stars). Finally, this bias is important in our discussion 
only if cluster cr's are affected differentially as a function 
of age. 

Binaries are a challenging solution to the current prob- 
lem because such a model requires careful coordination of 
the binary fraction as a function of age, given the smooth 
rise in Y* around ages of 1 Gyr that mimics evolution, 
and then the subsequent lack of observationally impor- 
tant binaries in older systems. Despite such an argument, 
there is precedent for stellar evolution highlighting cer- 
tain stellar population s, with different bin ary fractions, 
at different times (see iGieles et alJl2010bl ). That work, 
however, focused on a population that was 10 Myr old, 
where it is easier to highlight short-lived populations. At 
ages of 1 Gyr a wider mass range of stars contribute sig- 
nificantly to the luminosity. We believe such a scenario 
unlikely, but cannot yet reject it categorically (see §6 for 
how we expect to rule it out). 

5.1. Are Evolutionary Models the Problem? 

The unexpected behavior of T* with age could reflect 
a problem with our expectations. Previously unappre- 
ciated complications that affect the calculatio ns of T* 
continue to be highlighted (|Fan fc de Griis 2012|). Due to 
the nature of the data in Figure [5l such an approach must 
predict falling values of T* with age for stellar popula- 
tions with ages somewhere between 9.4 < log(age [yrs]) 
< 9.8. Such a prediction does not pose a problem for our 
sample of clusters, which does not probe these epochs, 
but it might for galaxies, which do. Studies of early- 
type galaxies, which are presumed to have had relatively 
simple star formation histories, have concluded that the 
evolution out to z ~ 1 is broadly consistent with the 
expectation of passive evol ution for an intrinsically old 
population (f or exa mples seelyan Dokkum fc Franxll 19961: 
i Kelson et al.lll997tlTreu et al.ll2005allbl: Ivan de Ven et ahl 
120031: Ivan Dokkum k van der MarelP2007fi . These re- 
sults depend on observations from which the investiga- 
tors infer that values of Y* are dropping with increasing 
redshift, or decreasing age, rather than rising. 

To determine whether galaxies can help discriminate 
against strong swings in Y* at intermediate ages, we use 
measurements of the ages and structural parameters of 
two independent sets of early type galaxies. First, we 
use the study of Sloan Di gital Sky Su r vey ( SPSS) lo- 
cal galaxies presented by iGraves et al.l (|2009l ). In this 
study, they binned early-type galaxies by their struc- 
tural parameters, a, fh, and Ih, a nd examined devia- 
tions from the Fundament al Plane (jDressler et alJlT987t 
IDiorgovski k Davlsl I1987T) relative to age and chemi- 
cal abundance as measured from the stacked spectra of 
galaxies in that bin using the Balmer absorption lines 
and Lick indexes. These are luminosity weighted quan- 
tities and subject to the usual caveats involved in pa- 
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rameterizing complex populations by a single age and 
metallicity. Nevertheless, we use our Equation [2] to 
calculate values ol and compare those values, as a 
lunction of age, to those of our stellar clusters. One 
important distinction is that unlike the clusters, these 
galaxies do contain dark matter and the exact propor- 
tion of dark matter within is unknown and likely to 
vary as a function of these structural parameters (for 
some examples from the long history of this topic see 
Babul fc Reeslll99l iGrahaml 12001 iMarinoni fc Hudson! 
2002t Ivan den Bosch et al.ll2003t ICappellari et al.l l2006t 
Zaritsky et al.ll2006HWolf et al.ll2010( ). The results of our 
calculation are pre sented in Figured Second, we use the 
data presented by Ivan Dokkum fc van der Marell (f2007h 
who ascribed deviations from the Fundamental Plane for 
samples of galaxies in galaxy clusters as a function of 
redshift to evolution in T*. For these galaxies, we do 
not have a directly measured age estimate, so instead 
we assume that these galaxies formed at high redshift 
(here we take the formation time to be log(age) = 10.11, 
so that an elliptical galaxy at z — has the same age 
as our old stellar clusters) and that they have evolved 
passively thereafter. The age corresponding to a partic- 
ular redshift is calculated using Ho = 70 km s _1 Mpc -1 , 
Sl m = 0.3, and Q\ = 0.7. We correct their values of 
Y», which are in the rest B-band, to V by using the 
calculated color evolution of a stellar population with 
an instantaneous burst of star formation at log(age) = 
10.11 using PEGASE (jFioc fc Rocca-Volmerange1ll997| ) 
and adopting default values for th e binary fractio n (0.5) , 
stellar wind prescriptions, and a iMiller fc Scald (|1979f ) 
IMF that extends between 0.1 and 120 M Q . Because 
the measured shifts from the z — Fundamental Plane 
do not provide us with an absolute calibration of T*, 
the authors tabulate values for AT* as a function of 
z, there is some freedom in normalization. We normal- 
ize the value s of T h to match those calculated for the 
iGraves et al.l (|2009l ) galaxies, al though an alterna t e nor - 
malization using the results of ICappellari et al.l (|2006f ) 
places the I van D okku m fc van der Marell (|2007| ) slightly 
below the IGraves et al.l (j2009T ) galaxies. In either case, 
the excellent agreement in the slope of the trend of Th 
with age is important because the galaxy ages in these 
two samples are measured so differently. 

These results (Figure [9j vertical axis changed to 
log^T*) for a reason that will become apparent below) 
represent the upper bounds of T* in early-type galaxies 
because a correction for the contribution of dark matter 
to Y/,, will drive T* downward. The exact magnitude 
of the correction is uncertain because it depends on the 
very stellar evolutionary models that we are attempting 
to test. We can posit that dark matter is responsible for 
the entire offset between this population and our old stel- 
lar clusters, if we hypothesize that the stars in early-type 
galaxies are drawn from the same parent population as 
those in clusters. To arrive at this agreement would re- 
quire an 80% dark matter mass fraction, and we show the 
resulting comparison also in Figure [9l Perfect agreement 
is unlikely because the galaxies are less likely than the 
clusters to be described well as the result of instantaneous 
burst of star formation and the dark matter fraction will 
vary among galaxies. 

In the same Figure, we also include two model tracks 
for the evolution of simple stellar populations. First, we 
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Fig. 9. — The evolution of T* as inferred from early type galax- 
ies and comparison to models and our cluster data. Our clus- 
ter data are plotted i n solid circles (black), the resu lts obtained 
using the results from van Dokkum & Conroy (2010) are plotted 
with open cir cles (red) when norm alized to match the results from 
the sample of Graves et al. ( 2009 ), which are plotted with squares 
(pink). We correct the van Dokkum & Conroy (2010) results to the 
V band using a simple stellar population model, but no correction 
for dark matter (as such these ar e maximal T*). We also show 
the van Dokkum & Conrov (2010) data assuming a dark matter 
fraction of 80%, so as to qualitatively match to our stellar clus- 
ter data, in open circles (blue). The likely dark matter fraction is 
likely to be somewhere between and 80%. The upper solid line 
represents the values of Tj, from a PEGASE model of a population 
with an instantaneous burst at t = and a Salpeter IMF. while 
the lower represents a model with a light- weighted Kroupa IMF. In 
the log-log space used here, the model tracks can be approximated 
as straight lines of slope —0.77. 

note that we have corrected the plotted values of T* for 
the stellar clusters for the effect o f internal dynamica l 
evolution using the results of the lAnders et al.l (|2009l ) 
models. We plot the results of PEGASE models using a 
Salpeter IMF, spanning from 0.1 to 120 Mq, with default 
stellar mass loss and binarity parameters, and metallic- 
ity matching the mean of the young clusters (—0.4). We 
have not renormalized these results in any way and they 
do an acceptable job of reproducing the trend seen in 
the younger clusters, for clusters of 8 <log(age [yrs]) < 
9.4. As we discussed previously, we ignore discrepancies 
at younger ages due to the possibility that these clus- 
ters are not relaxed. The slight systematic underesti- 
mation of T* could be addressed by extending the IMF 
to somewhat lower mass and perhaps adding some sub- 
stellar mass, which is set to zero in the current models. 
We also plot results obtained using a Kroupa IMF, with 
the same parameters as the Salpeter model except that 
the metallicity is set to match that of the older clusters 
(—1.8) and here we are forced to renormalize the results 
by removing half of the stellar mass to reach agreement 
with the properties of the older clusters. We therefore 
refer to the function shown as a light-weighted Kroupa 
IMF. 

The comparisons between clusters, galaxies and mod- 
els provide several interesting results. First, we confirm 
that the early-type galaxies require that T* decrease with 
lookback time rather than increase over the epoch where 
we lack cluster measurements. This result rules out mod- 
els with strong swings in T* at intermediate ages. Sec- 
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ond, the early- type galaxies lie between the extrapolation 
of the young cluster trend and the old clusters. This "dis- 
agreement" is is easily addressed because the measured 
values of are not equivalent to Y* due to dark mat- 
ter's contribution to the mass. Third, we are almost at 
the point with the galaxy sample where it overlaps in age 
with the intermediate age clusters. Such overlap would 
provide a valuable test of the models because even though 
the two types of systems appear at the same age, they are 
inherently different in that for the clusters we are plot- 
ting stars formed X years ago, while for the galaxies we 
are plotting stellar populations that presumably formed 
13 Gyr ago but that we are seeing as they were X years 
after their formation. Fourth, we find another reason for 
dismissing dramatic dynamical evolution of the old clus- 
ters as a solution to our problem because the values of 
for the galaxies already lie below the extrapolation 
of the young clusters, and their T* values will lie even 
lower due to the dark matter correction. 

6. THE INITIAL MASS FUNCTIONS 

We conclude that the solution to the apparent discrep- 
ancy illustrated by Figures |6] and |9] is most directly ad- 
dressed with the existence of at least two distinct stellar 
initial mass functions. 

Although our data are limited by sample size, they 
can begin to inform a discussion regarding the cause 
of IMF variations. For example, one plausible candi- 
date for the driver of those variations is metallicity, al- 
though in fact little variation in the IMF has been seen to 
date with metallicity and some theoretical under standing 
for that finding has emerg ed ()Krumholdl20rl . In our 
data, however, all of the old clusters are of low metal- 
licity ((Fe/H) 2) and NGC 416, which is the only 

cluster with log(age) < 10 that lies with the old clus- 
ters is also of relatively low metallicity (—1.44). NGC 
416 is a particularly interesting candidate in this con- 
text because it formed well after the initial episode of 
star forma tion in the SMC (being only 7 Gyr o ld). If we 
search the M cLaughlin fc van der Marell (f2005h data for 
other clusters that are of intermediate age 9 < log (age 
[yrs]) < 10 and low metallicity ((Fe/H) < —1), we find 
only three additional such clusters (KRON 3, NGC 339, 
and NGC 361). Unfortunately, the mean surface bright- 
ness within th of these clusters is at least 7 times fainter 
than that of NGC 416, which already required over 5000 
sec of exposure time. An alternative test is to find 
old clusters (log(age [yrs]) > 10) with high abundances 
((Fe/H) > —0.5). In this case we find 11 clusters (Liller 
1, NGC 5927, NGC 6440, NGC 6528, NGC 6553, NGC 
6624, Pal 8, Pal 10, Pal 11, Terzan 2, and Terzan 5). Of 
these, NGC 6440, 6528, 6553, and Pal 11, have the re- 
quired data, other than cr, and high mean surface bright- 
ness, and are therefore good candidates for subsequent 
study. 

To uncover any additional constraints on potential 
drivers of di stinct IMFs, we combine our data a nd that 
compiled by [McLaughlin fc van der Marell (|2005T ). In the 
log-log space of Figure GO the model tracks are roughly 
linear with slope —0.77. As such, it is straightforward to 
predict the value of T* for each cluster at any age. We 
choose to calculate Y* at 10 Gyr, Y^^o, for all of the clus- 
ters that are older than log (age [yrs]) > 8.3 as a way to 
directly compare them. The differences among values of 



Y*,io will reflect differences in the IMF, if (1) the effects 
of dynamical evolution have been correctly accounted for, 
(2) stellar evolution models do follow the adopted linear 
behavior (in log space), and (3) clusters along any partic- 
ular evolutionary track are of similar chemical abundance 
or chemical abundance does not grossly affect Y^io- We 
present the distribution of Y^io versus age and chemical 
abundance in Figure fTOl 
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Fig. 10. — Evidence for two initial mass functions. We plot 
Y^iOi the value of T* at 10 Gyr for each cluster obtained as 
described in the text, versus iron abundance and age. The clus- 
ters included are those from our sa mple (with errorbars) and the 
IMcLaughlin fc van der Mare l ( 2005) sample for log(age) > 8.3. For 
the McLaughlin & van der Marel (2005) clusters we adopt their 
calculated values of T* from the dynamical models, except we re- 
duce those esti mate s by a factor of 0.7 to match normalizations to 
our estimates ( jQTT) - For our clusters, we adopt the values of Y„ 
tha t include the correct ions for dynamical evolution obtained from 
the I Anders e t al. (2009) models. The shaded regions represent the 
single cluster rms scatter of each population, centered on the mean 
value of Y^io for each population. 

The distribution of T^io is striking. There are two 
well-separated populations of clusters, highlighted by the 
grey bands that represent the mean and single cluster 
standard deviation for each of the two populations. The 
means of each population are statistically distinct at 
~ 10cr. It is also evident that neither age nor (Fe/H) 
appear to be the sole arbiter of stellar population proper- 
ties. Unfortunately, neither the clusters (NGC 2257 and 
6535) at low (Fe/H) and high Y^io, which best overlap 
the older clusters, nor the low (Fe/H), low Y* 10 clus- 
ters (NGC 104, 6362, 6388, and 6441, excluding NGC 
6366 which has an anomalously low measurement of Y„) 
that overlap the young population are among those for 
which we obtained measurements of the internal kine- 
matics. As such, we cannot determine the confidence 
with which these are clusters that provide evidence for 
overlapping properties among these two populations. If 
their velocity dispersions are confirmed by our upcom- 
ing measurements, then these clusters would be key in 
unraveling the nature of the two populations. 

Direct measures of the mass function, obtained by 
counting stars, is evidently the most robust way to de- 
termine if there are indeed IMF variations among the 



14 



Zaritsky et al. 



clusters. Unfortunately, the Galactic clusters for which 
this has been done extremely well are all old clusters, 
and therefore lie predominantly in only one of our two 
populations. Nevertheless, these results can still shed 
some light on some of our results. First, for NGC 6366, 
which is the outlier with extremely low T^io in Figure 
[TU1 the direct measurement of the stellar mass function 
from Hubble Space Telescope imaging demonstrates that 
its mass function is indeed anomalous in th at it has far 
fewer low mass stars than other clusters (jPaust et al.l 
2009). This results at least qualitatively valida tes the low 
value shown in the Figure. IPaust et al.l (|2009D attribute 
the low number of low mass stars, or correspondingly our 
low value of T^io, to strong tidal stripping because this 
cluster lies in the Galactic bulge. Second, excluding NGC 
6366, the range of Y» among the group of low T* clusters 
is roughly 0.5 dex or a factor of 3. Using HST imaging 
of a set of Galactic clusters (w hich are all old an d so 
presumably fall into this group) . IPaust et all (|2010[ ) find 
that the stellar mass function slopes, for M < 0.8M Q , 
vary from —1.7 to —0.3, where the Salpeter slope is —2.35 
in this convention. This range of slopes corresponds to 
a factor of 3.6 change in the mass contained in stars be- 
tween 0.1 and 0.8 M Q , and so to a somewhat smaller 
change in the overall mass. We conclude that much of 
the apparent scatter in T^io could be due to these slope 
differences, which are presumably due to different inter- 
nal evolution. This finding is further confir med with a 
sample of SMC clusters (jGlatt et al.ll2011| ). Here the 
clusters span ages from 9.8< log (age [yrs]) < 10.2, so 
again they unfortunately do not reach the clusters in our 
high Y* population, but have slopes ranging from —1.4 
to —1.74, if one excludes two ambiguous determinations 
and a cluster that we discuss below. The two ambiguous 
results are for NGC 121 and NGC 416, two clusters that 
are also part of our sample, for which the slope is either 
around —1.4, which would be entirely consistent with the 
Galactic clusters, or —2.3, depending on whether the cen- 
tral region is excluded or replaced with higher resolution 
HST HRC data. In one case the slope decreases when 
the inner data are included, in the other it increases, so 
the uncertainty here is large. A third cluster, NGC 339, 
with a steep mass function slope is potentially an inter- 
esting target because that one might lie in the high T* 
population, although there is no measurement of a yet 
available. 

7. SUMMARY 

Using determinations of the velocity dispersions of 18 
stellar clusters within the Local Group that span 7 < log 
(age [yrs]) < 10.2, we trace the evolution of the stellar 
mass-to-light ratio of simple stellar populations. We find 
that the observed behavior does not match simple the- 
oretical expectations, and that the failure is such that 
global scalings of the masses, luminosities, or universal 
initial mass function cannot reconcile the two. The ef- 
fect of the internal dynamical evolution of the clusters on 
T* goes in the correct sense to alleviate the discrepancy, 
but fails to produce results that the sufficiently large. 
This is confirmed both by using models of the dynami- 
cal evolution (jAnders et al.ll2009l ) and direct measures of 
the s t ellar mass function in nearby clusters ()Paust et al.l 
2010: iLeigh et al.ll2012t) . Contamination by binaries is 



also discussed, but seems problematic due to fine tuning 
problems. 

The data can be explained if there are two populations 
of clusters, drawn from different initial mass functions. 
Specifically, we find that one IMF is primarily, but not 
exclusively, appropriate for older, metal poor clusters and 
the other for primarily, but not exclusively, for younger, 
metal rich clusters. The young (log(age [yrs])<9.5) clus- 
ters are well-described by a bottom-heavy IMF, such as 
a Salpeter IMF, while the older clusters are better de- 
scribed by a top-heavy IMF, such as a light-weighted 
Kroupa IMF, although neither of these specific forms is a 
unique solution. Of course, the sample is currently small 
and there are only four clusters in the key age range. 
We are in the process of enlarging the sample, although 
there are not many intermediate age clusters available 
for such study with the current technology. There are 
also a variety of effects, including the dynamical evolu- 
tion described above, whose magnitude is estimated us- 
ing models; therefore, significant departures from those 
models can also be used to explain the observations. 

Claims for variations in the initial mass function are 
usually treated with significant skepticism. This is not 
because we have a basic understanding of the initial mass 
function that predicts universality, but is rather due to 
a fondness for Occam's razor. In fact, certain previous 
claims of variations in the IMF ar e in the same s ense as 
the effect seen here. For example, iKroupal (|2001| ) found 
evidence for greater num bers of low m ass stars in younger 
Galactic populations and Dave (2008) inferred top heavy 
star formation at earlier times on the basis of galaxy evo- 
lution models. Whether these findings are quantitatively 
consistent with what we present is unclear, but certainly 
worth exploring to determine if a consistent picture can 
emerge. 

It is certainly the case that allowing for the possibility 
of varying IMFs tremendously complicates what one can 
constrain in extragalactic astronomy. Even if there are 
only two IMFs, the variation between those as a function 
of age, environment, star formation rate, or whatever else 
drives the dichotomy provides numerous new degrees of 
freedom to any simulation or model. Of course, some 
have already appealed to such effects to alleviate conflicts 
between models and data, but multiple IMFs not only 
enables one to address what were previously frustrating 
conflicts, but requires one to reassess what were previ- 
ously apparent successes. Even worse, without knowing 
the underlying physical cause for the two classes of IMFs 
claimed to be responsible for the observations presented 
here, we cannot claim that there are only two variants. 
Unfortunately, this is now an added complication that 
we must face or at least include as a source of systematic 
error in any study involving galaxy luminosities, color, 
or stellar mass determinations. 



DZ acknowledges financial support from NSF grant 
AST-0907771 and NASA ADAP NNX12AE27G and 
thanks the Max Planck Institute for Astronomy for its 
hospitality during the completion of this study. R. C. ac- 
knowledges support from NSF through CAREER award 
0847467. 



15 



REFERENCES 



Anders, P., Lamers, H. J.G.L.M., & Baumgardt, H. 2009, A&A, 
502, 817 

Aragon-Salamanca, A., Baugh, CM., & KaufTmann, G. 1998, 

MNRAS, 297, 427 
Babul, A., & Rees, M.J. 1992, MNRAS, 255, 346 
Bahcall, J. N., & Tremaine, S. 1981, ApJ, 244, 805 
Bastian, N., Covey, K.R., & Meyer, M.R. 2010, Ann. Rev. 

Astron. Astroph., 48, 339 
Baumgardt, H., & Makino, J. 2003, MNRAS, 340, 227 
Bell, E.F., Mcintosh, D.H., Katz, N., & Wienberg, M.D. 2003, 

ApJS, 149, 289 
Bell, E.F., et al. 2004, ApJ, 608, 752 
Bender, R. 1990, A&A, 229, 441 

Bernstein, R., Shcctman, S.A., Gunnels, S.M., Mochnacki, S., & 

Athey, A.E. 2003, SPIE, 4841, 1694 
Bruzual, G. & Chariot, S., 2003, MNRAS, 344, 1000 
Bundy, K. et al. 2006, ApJ, 651, 120 
Cappellari, M. et al. 2006, MNRAS, 366,1126 
Cappellari, M. et al. 2012, arXiv:1208:3523 

Chandar, R., Fall, M.S., & Whitmore, B.C. 2010, ApJ, 711, 1263 
Chabrier, G., 2003, PASP, 115, 763 

Colucci, J.E., Bernstein, R.A., Cameron, S.A., & McWilliam, A. 

2011, ApJ, 735, 55 

Colucci, J.E., Bernstein, R.A., Camecon, S.A., & McWilliam, A. 

2012, ApJ, 746, 29 

Conroy, C. & Gunn, J. E. & White, M., 2009,MNRAS, 699, 486 
Dabringhausen, J., Hilker, M., & Kroupa, P. 2008, MNRAS, 386, 
864 

Dabringhausen, J., Kroupa, P., Baumgardt, H. 2009, MNRAS, 
394, 1529 

Dalcanton, J. J., et al. 2012, ApJS, 198, 6 

Dave, R. 2008, MNRAS, 385, 147 

Daddi, E., et al. 2005, ApJ, 626, 680 

Djorgovski, S. & Davis, M., 1987, ApJ, 313, 59 

Dressier, A. & Lynden-Bell, D. & Burstein, D. & Davies, R. L. & 

Faber, S. M. & Terlevich, R. & Wegner, G., 1987, ApJ, 313, 42 
Dutton, A., Mendel, J.T., & Simard, L. 2012, MNRAS, 422, 33 
Fan, Z., and de Grijs, R. 2012, MNRAS, 424 2009 
Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Gieles, M. Baumgardt, H., Heggie, D.C., & Lamers, H.J.G.L.M. 

2010, MNRAS, 408, L16 
Gieles, M. Heggie, D.C., & Zhao, H.-S. 2011, MNRAS, 413, 2509 
Gieles, M., Sana, H., & Portegies Zwart, S.F. 2010, MNRAS, 402, 

1750 

Giersz, M., & Heggie, D.C. 1997, MNRAS, 286, 709 
Glatt, K. et al. 2011, AJ, 142, 36 

Goodwin, S. P., & Bastian, N. 2006, MNRAS, 373, 752 

Graves, G.J., Faber, S.M., & Schiavon, R.P 2009, ApJ, 698, 1590 

Graham, A. 2002, MNRAS, 334, 721 

Henon, M. 1965, Ann d'Astrophys. 28, 62 

Henon, M. 1969, A&A, 2, 151 

Henon, M. 1975, in IAU Symp. 69, (cd. Hayli, A.), Reidel, 

Dordrecht, p. 133 
Hills, J. G. 1980, ApJ, 235, 986 
Jordi, K. et al. 2009, AJ, 137, 4586 
Juneau, S. et al. 2005, ApJ, 619, 135 

Kaufmann, G., White, S. D. M., Heckman, T., Menard, B, 

Brinchmann, J., Chariot, S., Tremonti, C, & Brinkmann, J. 

2004, MNRAS, 353, 713 
Kelson, D., van Dokkum, P.G., Franx, M., Illingworth, G., & 

Fabricant, D. 1997, ApJ, 478, 13 
Kroupa, P., Tout, C. A., Gilmore, G., 1993, MNRAS, 262, 545 
Kroupa, P. 2001, MNRAS, 322, 231 
Kruijssen, J.M.D. 2008, A&A, 486, L21 
Kruijssen, J.M.D. & Mieske, 2009, A&A, 500, 785 
Krumholz, M.R. 2011, ApJ, 743, 110 
Kuijken, K. & Merrifield, M.R. 1993, MNRAS, 264, 712 
Lane, R.R., Kiss, L.L., Lewis, G.F., Ibata, R.A., Sibert, A., 

Bedding, T.R., & Szckcly, P. 2010, MNRAS, 401, 2521L 
Langer, N., & Maeder, A. 1995, A&A, 295, L685 
Lamers, H. J. G. L. M., Gieles, M., Bastian, N., Baumgardt, H., 

Kharchcnko, N. V., & Portegeis Zwart, S. 2005, A&A, 441, 

117L 



Lamers, H J. G. L. M., Gieles, M., Zwart, S.F.P. 2005, A&A, 429, 
173L 

Leigh, N., Umbreit, S., Sills, A., Knigge, C, de Marchi, G., 

Glebbeek, E., & Sarajedini, A. 2012, MNRAS, 422, L1592 
Maraston, C. & Daddi, E. & Renzini, A. & Cimatti, A. & 

Dickinson, M. & Papovich, C. & Pasquali, A. & Pirzkal, N., 

2006, ApJ, 652, 85 
Marinoni, C, & Hudson, M.J. 2002, ApJ, 569, 101 
McLaughlin, D. E. & van der Marel, R. P., 2005, ApJS, 161, 304 
McQuinn, K.B.W., Skillman, E.D., Dalcanton, J. J., Dolphin, 

A.E., Holtzman, J., Weisz, D., & Williams, B.F. 2011, ApJ, 

740, 48 

Melbourne, J., et al. 2012, ApJ, 748, 47 

Meylan, G., Dubath, P., & Mayor, M. 1991, IAU 148 (Kluwer: 

Dordrecht), p. 211 
Miller, G. E., & Scalo, J. M. 1979, ApJS, 41, 513 
Milone, A.P. et al. 2012, A&A, 540 16 

Minor, Q.E., Martinez, G., Bullock, J., Kaplinghat, M., & 

Trainor, R. 2010, ApJ, 721, 1142 
Narayanan, D., & Dave, R. 2012, MNRAS, 423, 3601 
Olszewski, E.W., Pryor, C. & Armandroff, T.E. 1996, AJ, 111, 

750 

Padmanabhan, N. et al. 2004, New Astronomy, 9, 329 
Paust, N.E.Q., et al. 2009, AJ, 137, 246 
Paust, N. et al. 2010, AJ, 139, 476 
Popescu, B, & Hanson, M. M. 2010a, ApJ, 713, 21 
Popescu, B, & Hanson, M. M. 2010b, ApJ, 724, 296 
Portegies Zwart, S.F., McMillan, S.L.W., & Gieles, M. 2010, 
ARAA, 48, 431 

Pryor, CP., Latham, D.W., & Hazen, M.L. 1988, AJ, 96,123 
Rafclski, M., & Zaritsky, D. 2005, AJ, 129, 2701 
Rejkuba, M., Dubath, P., Minniti, D., & Meylan, G. 2007, A&A, 
469, 147 

Rieke, G.H., Loken, K., Rieke, M.J., & Tamblyn, P. 1993, ApJ, 
412 99 

Rix, H.-W., & White, S.D.M. 1992, MNRAS, 254, 389 
Rudnick, G. et al. 2003, ApJ, 599, 847 
Salpcter, E.E. 1955, ApJ, 121, 161 

Sargent, W.L.W., Schechter, P.L., Bokensberg, A., & Shortridgc, 

K. 1977, ApJ, 212, 326 
Spitzcr, L. 1958, ApJ, 127, 17 
Spitzer, L., & Mathieu, R. D. 1980, ApJ, 241, 618 
Spitzer, L., Jr. & Shull, J. M. 1975, ApJ, 201, 773 
Thomas, D., Maraston, C, Bender, R., & Mendcs de Oliveira, C. 

2005, ApJ, 621, 673 
Tremonti, C.A. et al. 2004, ApJ, 613, 898 

Treu, T. Ellis, R.S., Liao, T.X., van Dokkum, P.G. 2005, ApJ, 622, 
5 

Treu, T. et al. 2005, ApJ, 633, 174 

Treu, T., Auger, M.W., Koopmans, L.V.E., Gavazzi, R., 
Marshall, P.J., & Bolton, A.S. 2010, ApJ, 709, 1195 

van de Ven, C, van dokkum, P.G., & Franx, M. 2003, MNRAS, 
344, 924 

van den Bosch, F. C, Yang, X., Mo, H.J. 2003, MNRAS, 340,771 
van Dokkum, P., 2008, ApJ, 674, 29 

van Dokkum, P.G., & Conroy, C. 2010, Nature, 468, 940 

van Dokkum, P.G. & Franx, M. 1996, MNRAS, 281, 985 

van Dokkum, P. C, & van der Marel, R. P. 2007, ApJ, 655, 30 

Walker, M.G., Mateo, M., Olszewski, E.W., Penarrubia, J., Wyn 

Evans, N., & Gilmore, G. ApJ, 704, 1274 
Waters, CZ, Zepf, S.E., Lauer, T.R., Baltz, E.A., & Silk, J. 2006, 

ApJ, 650, 885 
Whilcy, I.M. et al. 2008, MNRAS, 387, 1253 

Whitmore, B.C., Chandar, R., & Fall, S.M. 2007, AJ, 133, 1067 
Winsall, M.L. & Freeman, K.C 1993, A&A, 268, 443 
Wilson, CP. 1975, AJ, 80, 175 

Wolf, J., Martinez, G.D., Bullock, J.S., Kaplinghat, M., Geha, M., 
Munoz, R.R., Simon, J.D., & Avedo, F.F. MNRAS, 406, 1220 

Zaritsky, D. & Gonzalez, A. H. & Zabludoff, A. I., 2006, ApJ, 
638, 725 

Zaritsky, D. & Zabludoff, A. I. & Gonzalez, A. H., 2008, ApJ, 
682, 68 

Zaritsky, D. & Zabludoff, A. I. & Gonzalez, A. H., 2011, ApJ, 
727, 116 



